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Constraining a fourth generation of quarks: 
non-perturbative Higgs boson mass bounds 
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We present a non-perturbative determination of the upper and lower Higgs boson mass bounds 
with a heavy fourth generation of quarks from numerical lattice computations in a chirally symmetric 
Higgs- Yukawa model. We find that the upper bound only moderately rises with the quark mass 
while the lower bound increases significantly, providing additional constraints on the existence of a 
straight-forward fourth quark generation. We examine the stability of the lower bound under the 
addition of a higher dimensional operator to the scalar field potential using perturbation theory, 
demonstrating that it is not significantly altered for small values of the coupling of this operator. 
For a Higgs boson mass of ~ 125GeV we find that the maximum value of the fourth generation 
quark mass is ~ 300GeV, which is already in confiict with bounds from direct searches. 
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INTRODUCTION 

A heavy fourth generation of quarks is an attrac- 
tive and simple extension (denoted SM4) of the three- 
generation Standard Model (SMS) However, heavy 
fermion effects are expected to significantly contribute 
to the properties of the Higgs boson, leading to measur- 
able deviations in Higgs production cross sections and 
branching ratios. The recent discovery of a scalar bo- 
son seemingly consistent with the Standard Model expec- 
tation therefore casts serious doubt on straight-forward 
fourth generation scenarios [s^ . 

Another property of the Higgs boson which is sensitive 
to a possible fourth generation is its mass. Invoking ar- 
guments from perturbation theory, the Higgs boson mass 
is bounded from above by the triviality bound, which re- 
flects the Gaussian nature of the UV fixed point, and 
from below by the vacuum instability bound, which en- 
sures that the theory has a stable vacuum state. 

In perturbation theory, the lower bound can be ob- 
tained by examining the effective potential and demand- 
ing that it remains bounded from below. As the fermion 
fields contribute negatively to the effective potential, they 
have a destabilizing effect which leads (by demanding 
the stability of the theory) to a lower Higgs boson mass 
bound. However it is expected that the perturbative ex- 
pansion breaks down for Yukawa couplings near or less 
than the perturbative unitarity bound which allows 
maximal fermion masses of roughly 
TO/ - 500-600GeV 0,0. 

Due to these considerations, study of heavy fourth gen- 
eration extensions to the Standard Model necessitates 
non-perturbative methods. To this end, we employ lat- 
tice field theory techniques which allow us to compute 



lower and upper Higgs boson mass bounds as well as reso- 
nance parameters of the Higgs boson non-perturbatively. 
Such a strategy has already been applied in Refs. .7.~10| 
for the non-perturbative determination of the upper and 
lower Higgs boson mass bounds and the Higgs boson res- 
onance parameters in the SMS. First results for a fourth 
quark generation have been presented in Ref. 11| . Inves- 



tigations of a bulk-phase transition at very large values 
of the Yukawa coupling have been initiated in Ref. [l^l 
while first studies of the system at non-zero temperatures 
can be found in Ref. [l3|. For a summary of recent work 
on this model see Ref. jl4{ . 

The lattice calculations are possible due to a lattice dis- 
cretization of the fermion action which respects an exact 
chiral symmetry at finite lattice spacing a 15[, thus en- 
suring that the chiral character of the Higgs-fermion cou- 
plings is respected by the lattice regulator in a conceptu- 
ally clean manner. This advance has triggered a number 
of lattice investigations of Higgs- Yukawa like models 0- 

In this letter we report on a systematic investigation of 
the lower and upper Higgs boson mass bounds for quark 
masses ranging from the Standard Model top quark mass 
of TO/ — nit « 175GeV up to masses of to/ « 700GeV, 
which is above the perturbative unitarity bound. Our 
calculations are performed at a fixed lattice cutoff of A = 
i « 1.5TeV which ensures that the fermion and Higgs 
boson masses are sufficiently far from the cutoff scale. 
We also rely heavily on the techniques and simulation 
strategies of Refs. [THlot. 

Our results indicate that with increasing quark mass 
there is a substantial upward shift of the lower Higgs 
boson mass bound leading to severe constraints on a 
straight-forward fourth quark generation in light of the 
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recently discovered scalar boson when interpreted as the 
Standard Model Higgs boson. In order to examine the 
stability of our results, we analyze the effect of including 
a higher dimensional operator in the scalar field poten- 
tial using perturbation theory. For small values of the 
coupling of this term we confirm that the lower bound is 
not significantly altered. 



THE MODEL 

Here we briefly review the definition of our model. For 
a more complete treatment, see Ref. Q. We employ a lat- 
tice discretization of the Standard Model Higgs-fermion 
sector which exactly respects a chiral symmetry at finite 
lattice spacing. As the dynamics of the complex scalar 
(Higgs) doublet are expected to be dominated by its in- 
teractions with the heaviest fermions, we consider a single 
degenerate quark doublet only. By the same reasoning, 
as well as for computational simplicity, we also neglect 
gauge fields. 

One of the defining features of the Standard Model is 
the chiral structure of the scalar-fermion interactions. In 
order to reproduce this structure in a lattice discretiza- 
tion it is crucial to maintain a chiral symmetry at finite 
lattice spacing. This has been a long-standing obstacle 
to the lattice regularization of Higgs- Yukawa models and 
was finally overcome by employing the Neuberger 'Over- 
lap' 15|, |22|, |23| discretization of the fermion action. 

Following the proposition in Ref. we can therefore 
construct a lattice Higgs- Yukawa model with a global 
SU{2)l X U{1)y symmetry. Specifically, the fields in- 
cluded in our model are a scalar doublet (f and two 
fermions, the left-handed components of which are paired 
into an SU{2) doublet. The lattice action can thus be 
written as 



S^Sf + S^, Sf ^'^■ipxMxyi'y, (1) 
xy 

+A^(c<fx^l)^ 

X 

Mxy = 2?i7'l2x2 + y{P+4P+ + P-^xP-)5xy, 

where is a doublet of four-component spinor fields, 
7)'^°^'^ is the free Overlap Dirac operator with a Wilson 
kernel, e IR'', and 4)^; = ^^t^, with = (1,-1^). 
The left- and right-handed projection operators P± and 
the modified projectors P± are given by 



P± = ^7^, P± = ^7^, 
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(2) 



The action introduced above obeys an exact global 
SU(2)i, X U{1)y lattice chiral symmetry. For € SU(2) 
and e [0, 27r] the action is invariant under the transfor- 
mation 



tP ^ UyP+'tlj + UYnLP^i;, 



(3) 



with Uy = cxp{i6Y), where Y labels the representation 
of the global hypercharge symmetry group C/(l)y. It 
should be noted that in the continuum limit the (global) 
continuum SU(2) ^ x U(1)y chiral symmetry is recovered. 

This formulation enables a numerical study of the limit 
A = oo on the lattice by simply enforcing the constraint 
^^^x = 1, Va;. Also, we can relate the parameters 
and fields appearing in Eq. [T] to those appearing in the 
standard scalar complex doublet continuum Lagrangian 
(/: = |9^^|2 + im2|^|2 + A|^|4)by 
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COMPUTATIONAL STRATEGY 

The cutoff in our model is provided by the inverse lat- 
tice spacing (A = 1/a) and we set the physical value of 
A using the phemonenological Higgs field vacuum expec- 
tation value (wfi), 



246 GeV = 



^ ■ a' 



(4) 



where Zq denotes the Goldstone boson field renormaliza- 
tion constant and v the bare scalar field vacuum expec- 
tation value. We target a mass range for the degenerate 
fermion doublet of 175GeV ^ to/ < 700GeV, while fix- 
ing the cutoff at A w 1.5TeV. Although the masses of 
the fermion doublet are degenerate in this work, we plan 
to assess the effect of a mass splitting in the near future. 
At m/ = TOt = 175GeV the splitting mf, — mt has been 
taken into account and found to have a small effect on 
the lower Higgs boson mass bound which, moreover, can 
be taken into account by the effective potential evaluated 
in lattice perturbation theory [Tj. 

We now briefly discuss the simulation algorithm. More 
details can be found in Ref. Q. The Monte-Carlo simu- 
lations are carried out using a variant [2?] of the Hybrid 
Monte Carlo (HMC) algorithm ^] in which the integra- 
tion over the Grassman flelds is done analytically at the 
expense of introducing the determinant of the fermion ac- 
tion in the updating. Therefore, our remaining degrees 
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of freedom in the Monte Carlo integration are the scalar 
fields. 

Due to the absence of gauge fields, the Overlap op- 
erator of Eq. ([ij can be constructed exactly in momen- 
tum space. However as the Yukawa coupling is local in 
position space, Fast-Fourier- Transforms (FFT's) are em- 
ployed to efficiently apply this term. Finally, Fourier 
acceleration [26[ is used to propagate the low momen- 
tum modes using coarser time-steps along the HMC evo- 
lution, which effectively reduces the autocorrelation be- 
tween successive scalar field configurations. 

It has been demonstrated Q that the Higgs boson mass 
is a monotonically increasing function of the bare quartic 
coupling A. Therefore, the lower bound for the Higgs 
boson mass at fixed cutoff and m/ is obtained at A = 0, 
while the upper bound is obtained at A = oo Q. 

Since we work in a finite volume with no external sym- 
metry breaking source, the naively defined vacuum ex- 
pectation value is zero in an ensemble average. We there- 
fore resort to a special techni que, pioneered in Ref. 
and employed in Refs. 
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of rotating a given scalar 
field configuration to a preferred direction. It can be 
shown that in infinite volume this leads to the same vac- 
uum expectation value as obtained in the standard proce- 
dure involving the limit of vanishing external source [30| . 

The determination of the Goldstone boson renormal- 
ization constant Zg, the Higgs boson mass uih, and 
the fermion mass m/ have been discussed in detail in 
Ref. and are briefly reviewed here. The renormal- 
ization constant Zq is computed from the slope of the 
inverse Goldstone boson propagator at vanishing Eu- 
clidean four-momentum transfer, in the standard on-shell 
scheme. Operationally, this constant is obtained from fits 
to the propagator at small momenta. 

Due to the existence of massless Goldstone modes in 
the spontaneously broken phase of our theory, the Higgs 
boson is unstable and can decay to final states contain- 
ing an even number of Goldstone bosons. We employ 
two definitions of the Higgs boson mass which both ig- 
nore this finite decay width. The Higgs boson propaga- 
tor mass is derived from fits to the momentum space 
Higgs propagator defined as 



hp = 



1 



while the correlator mass 



•'H 



is derived from fits to the 



Higgs temporal correlation function. Although finite de- 
cay width corrections to these formulae are quite differ- 
ent, the mass extracted from both of these procedures 
typically differs by less than 10%, lending credence to 
our approximation of a stable Higgs boson. Further- 
more, a rigorous study of the Higgs boson decay width 
at non-vanishing external source has been performed at 
175GeV in Ref. ^ which obtained a narrow 
decay width for all values of the bare quartic coupling, 



further supporting the validity of the stable Higgs boson 
approximation. For this work we quote as the central 
value as it is typically the lowest estimate of ttih ■ 

Finally, we compute the quark mass from the exponen- 
tial decay of the temporal correlation function C/(<) at 
large Euclidean time separations t, defined as 



LI _ 



(5) 



where the left- and right-handed spinors are defined using 
the projection operators in Eq. [2l 



NUMERICAL RESULTS 

As mentioned above, massless Goldstone excitations 
are present in the spontaneously broken phase of our 
model. Since we work at zero external symmetry- 
breaking source, finite size effects are not exponentially 
suppressed but rather algebraic in nature and propor- 
tional to inverse even powers of the linear extent of the 
lattice, i.e. C(l/L^) at leading order (isl-lsot. These fi- 
nite size effects can be quite substantial and an infinite 
volume extrapolation of our results for the renormalized 
vacuum expectation value and all masses is required. 

Finite volume data for the Higgs boson propagator 
mass (m|^) are shown in Fig. [T] as an example of our 
infinite volume extrapolations. The lattice data are plot- 
ted versus l/i^ and extrapolated to the infinite volume 
limit by means of a linear fit ansatz according to the 
aforementioned leading order behaviour. Due to the ob- 
served curvature arising from the non-leading finite vol- 
ume corrections, only those volumes with Ls > 16 have 
been included. One also finds that the infinite volume ex- 
trapolation can be reliably performed at ranges of lattice 



volumes from 12"^ x 32 to x 32 if subleading corrections 
are considered. 

The results of the Higgs boson masses for the lower 
and upper Higgs boson mass bounds as a function of the 
quark mass at a fixed cutoff of A = 1.5TeV are finally pre- 
sented in Fig. [21 All data for the upper bound have been 
extrapolated to infinite volume, as have the lower Higgs 
boson mass bounds at m/ « 160, 192, 420GeV Only lower 
Higgs boson mass bounds at to/ « 305, 500, 685GeV are 
taken on our presently largest lattice of size 24'^ x 32. 
However, the finite size corrections of these points will 
only provide a small change and the behaviour of the 
lower Higgs boson mass bound as a function of the 
fermion mass will not be significantly altered. 

We find that the upper Higgs boson mass bound is only 
moderately shifted by about 20% when the fermion mass 
is increased from 175GeV to 700GeV. On the contrary, 
the lower Higgs boson mass bound increases drastically 
with increasing quark mass. Such a significant shift of 
the lower bound has already been observed in f^, for a 
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FIG. 1: Infinite volume extrapolation of the Higgs boson mass 
extracted from fits to the momentum space Higgs propagator 
(m|J). The data are for the Higgs boson mass relevant for the 
upper bound (A = oo) and a range of bare Yukawa couplings 
corresponding to fermion masses from the physical top quark 
mass to ~ 700GeV. The infinite volume extrapolation is per- 
formed by fitting the data to a linear function in 1/ Ll taking 
only data with Ls > 16 into account. 
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FIG. 2: The lower and upper Higgs boson mass bounds 
as a function of the fermion mass nif for a fixed cutoff of 
A = 1.5TeV. The solid line represents a perturbative lattice 
effective potential calculation for the lower Higgs boson mass 
bound. 



The vacuum expectation value and lower Higgs boson 
mass bound are obtained from this effective potential in 
the usual way, namely 



(6) 



quark mass of m/ w 700GeV. We also show in Fig. [2] a 
lattice perturbative computation of the lower Higgs bo- 
son mass bound using the effective potential. The result 
of this calculation, which is discussed in the next sec- 
tion, is represented by the solid line in Fig. [2j Even up 
to fermion masses of my « 700GeV the lattice pertur- 
bative calculation qualitatively describes the simulation 
data rather well. This observation allows us to test the 
stability of the lower bound against the addition of a 
higher dimensional operator within the framework of the 
perturbative lattice effective potential. Moreover, direct 
Monte Carlo simulations with such a term were found 
to be well-described by perturbation theory Q within a 
small range of couplings of this higher dimensional oper- 
ator. 



EFFECTIVE POTENTIAL 

In addition to non-perturbative numerical data, it is in- 
structive to calculate the lower Higgs boson mass bound 
perturbatively using the effective potential. To this end 
we follow Ref. Q and calculate the effective potential to 
1-loop in the large- iV^: expansion. 

It should be noted that these calculations are per- 
formed with the same lattice regularization used in the 
simulations maintaining finite spatial and temporal ex- 
tents. Thus, the loop corrections are calculated numeri- 
cally, by explicitly performing sums over the lattice mo- 
menta. An infinite volume extrapolation of these pertur- 
bative results is also performed. 



where 4> denotes the average value of the scalar field. The 
results of such a determination are shown in Fig. [5J Al- 
though quantitatively different, the numerical data agree 
with the qualitative behaviour of the perturbative result 
even for rather large fermion masses. 

Based on this agreement, we estimate the effect of 
an additional higher-dimensional operator in the scalar 
field potential using the perturbative expansion discussed 
above. To this end we add to Eq. ^ the term Ag^^. It 
has been demonstrated Q that perturbation theory suc- 
cessfully describes non-perturbative numerical results for 
the range Xq = [0,0.01]. Therefore, at this preliminary 
stage we examine couplings in this range only. Further- 
more, in the presence of such a term we must modify the 
stability criterion of the effective potential. 

Due to the triviality of the theory, the cutoff A cannot 
be completely removed while maintaining non-zero val- 
ues of the renormalized quartic and Yukawa couplings. 
However, the existence of a scaling regime suggests that 
predictions of the model are affected only mildly by small 
cutoff effects provided that the cutoff is larger than the 
relevant scales in the theory, namely VR,mf, mg. It has 
been determined in the pure cj)'^ theory 3l| that a suit- 
able definition of the scaling regime is the situation where 
m/A < 0.5, where m = vii,mH- We adopt here this ob- 
servation also for m/, keeping mf/A < 0.5. 

It is in this scaling regime only that results which are 
universal up to small cutoff effects can be obtained from 
our model. Therefore, the vacuum stability criterion 
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should ensure that a stable vacuum exists throughout 
the scaling regime. Indeed, the running of the renormal- 
ized quartic and Yukawa couplings (and thus vacuum 
stability considerations) have been demonstrated to be 
severely regularization dependent outside of the scaling 
regime defined above [21|. Therefore, we choose as our 
stability criterion 



< 0.5. 



(7) 



In the Ag = case, this results in the well-known sta- 
bility criterion A > 0. We have examined the effect of a 
finite Ag in the range Ag = [0.0, 0.01] which results in a 
less than 1% difference in the lower bound. However for 
larger values of Ae ~ 0.1, deviations as large as 15% have 
been observed in the lower Higgs boson mass bound. As 
we are unsure of the applicability of perturbation theory 
in this region, further non-perturbative numerical simu- 
lations must be performed to properly assess the effect of 
this higher dimensional operator for larger values of Ag. 



OUTLOOK AND CONCLUSIONS 

In this letter we have performed a non-perturbative 
lattice investigation of the lower and upper Higgs boson 
mass bounds in an extension of the Standard Model with 
a heavy fourth quark generation. Since the heavy quark 
masses lead to large values of the Yukawa coupling, such 
a non-perturbative computation is necessary in order to 
have reliable results for the Higgs boson mass bounds at 
large m/. We include only the dominant interactions, 
which are expected to be the Higgs- Yukawa interactions 
of the heavy fermions. In particular, we neglect all gauge 
fields. 

Our chirally invariant lattice regularization of the 
Higgs- Yukawa sector of the Standard Model (as proposed 
in Ref. [l^]) obeys a global SU (2)l x U (1)y symmetry at 
finite lattice spacing, providing a formulation that pre- 
serves the chiral nature of the Higgs-fcrmion couplings. 
In this setup, for a fixed cutoff of A = 1.5TeV, we 
have studied a range of quark masses 175GeV ^ mf < 
700GeV. We find that the upper Higgs boson mass 
bound is only moderately increased by about 20% for 
larger quark masses. The lower Higgs boson mass bound 
however, changes quite substantially as is summarized in 
Fig. [2] and assumes for e.g. m/ « 500GeV a value of 
771^* « 500GeV. This puts severe constraints on the ex- 
istence of a fourth generation of quarks, given the current 
status of direct searches [H, |33| . 

Confronting the non-perturbatively computed lower 
Higgs boson mass bounds with a lattice perturbative 
calculation of the effective potential, we found that the 
effective potential describes the simulation data rather 
well, even for m/ w 700GeV. Based on this fact, as well 



as previous numerical work, we also used the perturba- 
tive effective potential to test the stability of the lower 
Higgs boson mass bound against addition of a higher di- 
mensional operator for a small range of couplings. We 
found that the lower Higgs boson mass bound is quite 
insensitive against such an addition, although more non- 
perturbative numerical simulations are required to exam- 
ine the effect of the higher dimension operator for larger 
couplings. Nonetheless, this preliminary work provides 
some confidence that the results for the lower bound are 
robust. Furthermore, although we work with a mass de- 
generate quark doublet in this work, the dependence of 
the lower bound on the mass splitting was found to be 
small at m/ = 175GeV, but will investigated at larger 
quark masses in the near future. 
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